Applied Probability Trust (19 October 2011) 



LOSS RATE FOR A GENERAL LEVY PROCESS WITH DOWN- 
WARD PERIODIC BARRIER. 

Z. PALMOWSKI,* University of Wroclaw 



1 



P. SWIATEK,** University of Wroclaw 



Abstract 

In this paper we consider a general Levy process X reflected at downward 
periodic barrier At and constant upper barrier K giving a process Vt^ = Xt-\- 
Lf — . We flnd the expression for a loss rate deflned by l^ = ELf" and 
identify its asymptotics as — >■ oo when X has light-tailed jumps and EXi < 0. 
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1. Introduction 

In this paper we consider a general Levy process X reflected at downward periodic 
barrier At — ip{t + U) (for a periodic nonnegative function ip{t) with period length 
s and U having uniform distribution on [0,s]) and at the constant upper barrier K, 
giving a process 

Vt"" ^ Xt + Lf ^ (1) 

being solution of respective Skorohod problem on each period where reflection is within 
bounded and convex set (see [51[T3). In above we assume that (p{t) G [0, a] for some 
a < K. Process X is defined on the filtered space (r2,-F, {J-t}t>o,V) with the natural 
filtration that satisfies the usual assumptions of right continuity and completion. From 
now on we will also assume that jump measure i/ of X is nonlattice. 
In this paper we find the expression for a loss rate defined by 

= ELf , (2) 

where E denotes the expectation when reflected process is stationary with stationary 
measure ttk (that is E(-) = Jq°°E[-|Vq'^ = x] nfcidx)), and prove that ~ De~^^ as 
K oo, where 7 solves ^(7) = for a Laplace exponent K(a) = logEexp{aX(l)} 
(which is well-deflned in some set Q) when X has light-tailed jumps and EXi < 0. 

The motivation comes from various queueing and telecommunication models (see 
[3 [31 in ini HI Hn [H]). Applications, where such a reflected Levy process considered in 
this paper is natural, are models where additionally to the input and output mechanism 
modeled by a Levy process there is a constant input given by a downward barrier At- 
This additional input is not available on liquid basis, but can only be used after some 
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maturity date has been reached. We choose that this time lag is fixed and equals 
the length s of the period of ip (for exponential time delay see [S]). For example, in 
view of internet networking applications, one considers the combined behavior of two 
services (e.g. streaming video and some other data). The first input behaves like a 
Levy process. The other input grows deterministically and can be served only at some 
fixed time s. The combined workload now behaves like a Levy process reflected at a 
lower barrier At. 

Fluid models with finite buffers are useful to model systems where losses are of 
crucial importance, as in inventory theory and telecommunications. Indeed, during 
the recent years real-time applications such as video streaming, interactive games 
have become increasingly popular among users. These applications are generally delay 
sensitive and require some preferential treatment in order to satisfy a desired level of 
Quality of Service. Traditionally finite capacity buffer mechanisms have been employed 
in the network routers, in which arriving packets are dropped when workload reaches 
its maximum capacity. In this paper we analyze the intensity of packets loss given by 
so-called loss rate given in 

There has been a great deal of work on overfiow probabilities in various fluid and 
queueing models but there have been relatively few studies on the loss rate in finite 
buffer systems. When jumps of the Levy process are heavy-tailed that one might have 
a hope to find relationship between these two notions (see [3j [6l [T^ for more classical 
models). Still in our model we focus on light-tailed case and therefore we choose Kella- 
Whitt martingale approach [S] . In fact we follow the ideas included in seminal paper of 
Asmussen and Pihlsgard [2], where both barriers are constant (see also [13] for matrix 
analytic method). This case corresponds to the assumption that A{t) = 0. Denoting 
the loss rate by analyzed in [5] we of course immediately get bounds: 



from which and [U Th. 4.1] it follows e.g. that \ogl^ — —7. In this paper we focus 
on more precise exact asymptotics. 

The paper is organized as follows. In Section [2] we give preliminary results and in 
Section [3] the main results with proofs. 



Assume from now on that (p G C^(int Jk) is invertible on some disjoint intervals Jk 
satisfying U)!^]^ Jfc = [0, s] with (p'{x) ^ for x e int Jk- 

Lemma 2.1. Process At — ip(t + U) has invariant measure: 



(3) 



2. Preliminaries 




where hk is a inverse of ip on int Jk ■ 



Proof. It is sufficient to check that for t £ [0, s): 



F{At <x)= F{ip{t + U) <x) = - [ 



s Jo 



'^{v{t+u)<x)du = - l(;^(„)<2.)du , (4) 
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where the last equahty is a consequence of periodicity of (p. The second part of 
the theorem follows from a straightforward arguments concerning a distribution of 
a piecewise strictly monotone function of r.v. 

□ 

Example 2.1. The most interesting case for applications is a saw-like lower boundary 
modeling constant intensity input (with rate 1 for simplicity): 

(p{t) = t mod a. (5) 

with < a < K . In this case n = 1, Ji = [0,a], s = a and £,{dy) = ^l(j,g[o.a])- 

Example 2.2. More complex situation will appear when ip is composed from few lines 
with different slopes: 

t for [0,1), 

V[t)^{ \-2{t-\) forte [1, 1), (6) 



3(t-|) forie[|,|). 



In this case rt = 3, Ji = [0,1], J2 = [1, f ), ^3 = [§,§), ^ ~ \ ^"^^ ^{dy) = 

i (i¥'(Ji)(y) + \^v{j2){y) + \\(Jz){y)) dy- 

Following arguments of Asmussen jX* P- 393-394] we have the following representa- 
tion of stationary distribution of . 

Lemma 2.2. The stationary distribution of the two-sided reflected Levy process 
is given by 



nv^ >x)= I > ' p{Xr > A-^"' + x)pu{z) adz), 



7r|^,(x) 



where t = inf{t > Q : Xt ^ [x ~ K. Af ^' ' -\- x)}, and A^ ^' = —(p(hk{z) — t) for 

z e ^(j,) andp,(z) = nu e j,\^{u) ^ z) = ■ 

Lemma 2.3. //E|Xi| < 00 then ELf < 00 and EL^ < 00 for each t>0. 

Proof Note that by ^ we have ELf < 00 if ELf < 00. Condition ELf < 00 
follows from the Wald identity applied to the random walk with increments being 
corrections of X between consecutive visits of the downward barrier (for details see 

HI). 

□ 

Further we need a little modification of Levy exponent «(«). We will treat large 
and small jumps separately. Let L be a constant that satisfies L > max(isr, 1). Then 
K{a) can be rewritten as follows: 

Oha+^^+j [e"^ - 1 -axl(|^|<i)]i/(da;), a G 9, (7) 
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where 6^ ^ 6 + xv{dx) + xv{dx). 

For any process Y we will denote by {Y^} its continuous part and the jumps by 

We split Aif^ into two parts, ALf and ALf , corresponding to AXs £ [0,L] 
and AXs G {L,oo), respectively, and ALf into AL^ and AL^, corresponding to 

AXs e [-L,0] and AAT^ e (-oo,-i), respectively. Let [f = EX;o<s<i4^f : = 

EEo<.<i ALf, 1/ = E Eo<.<i AL^, = E Eo<s<i Aif . Then /f = if +lf and 

^/ = if + Finally, let = - /f and 1^ = 1^- if with = ELf. 

Theorem 2.1. For a £ 6, 

Jo Jo 

+ ^ e»As (1 _ e--AL^ ) _ ^^<^K^K,c ^ ^ (1 _ e"^^") (8) 

0<s<t 0<s<t 

is a zero-mean martingale. 

Proof. It is well-known that for an J^t-adapted process = dY^ + X]o<s<t A^s 
of locally bounded variation, process 



Kt = K{a) f e'^^^ds + e"^ - e"^' +a f e"^=dr/ + ^ e"^=(l 



is a local martingale whenever a e 0, where = x + + Ft . Taking Yt — Lf — 
and using Lemma [231 to prove that Y has locally bounded variation, we get that 



Jo 

0<s<t 

-a f e^^'" dLf'^ + e"^"(l 



e 

0<s<t 



is a local martingale. Mt equals dH) since = K just after a jump of Lf, and 
= As just after a jump of Lf. To prove that {Mt} is a true martingale it sufficient 
to show that Esupo<s<t Ms < oo. This follows from the following conditions: < K, 

ELf•^ ELf''^ < cx) andEX;o<s<Jl-e"'^^^l < '^Y.o<s<t\'^ ~ e-'^^^'^l < °° (see 
also the proof of ^ Prop. 3.lJ). 

□ 
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Corollary 2.1. Let a G 0. Then satisfies the following equation: 
a(l - e"-^)Z^ = -K(a)Ee"^"'' + aEXi - ae"^l^ + off 

0<s<l 0<s<l 



0<s<l 0<s<l 

,2ti 



f AsdLf'" - a^E ^ A.ALf + o{a^) 

•'O 0<s<l 

Proof. If we take t = 1 in Mt and use the stationarity of V/^ we get 

= K(a)Ee"''«' + aE f e^-^^dif + E ^ e"^= (1 - e'^'^-^i') (10) 

0<s<l 



(9) 



ae 

0<s<l 



Moreover, we have 

^ (1 - e"^^^) = E (1 - + E (1 - ^"^'^)' (11) 



0<s<l 0<s<l 0<s<l 



0<s<l 0<s<l 0<s<l 

Applying expansion: 

= 1 + ax + + ^^e«"^ ^€[-1,1] (13) 

to the first parts of the right-hand sides of (ITT|) and and expansion e" = 1+q;-|-o(q;) 
to aE e°'^^dLf''^ completes the proof. 

□ 



We will need also the following observation. 
Lemma 2.4. We have: 



where EAq = J° y^(dy). 
Proof. Note that 



E / AsdLf = E(^o)^'^, 
Jo 



ELf = sEL^ = sl^ 
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which is consequence of the fact that Lf has independent and stationary increments 
under invariant starting position of A. The proof foUows now from the following: 



[ E [ A'JLf^idz) 
Jo 



adz) 



All' 



sl'^dA' 



adz) 



r lAll^^l^isAlll- f Aids) 
Jo L Jo 

pa pi pi pa 

/ / Atds^dz) ^l^ A 
Jo Jo Jo Jo 



adz) 

'■i{dz)ds = /^EAq. 



□ 



Now using Lemma we can rewrite © as follows. 
Lemma 2.5. As a ^0 we have 

a(l - e"-^ + aE{Ao))l^ = -K{a)Ee°'^'^ + aEX^ - ae°''^l^ + od^ 

+ ^E (Aif + E {AL^ 

0<s<l 0<s<l 



-e"^E J2 (l-e"^^=)-E ^ e"'^=(l 

0<s<l 0<s<l 

WE(Ao)EXi + a^E AsALf + 0(0^). (14) 

0<s<l 

3. Main results 

The first main result gives representation of in terms of basic characteristics of 
the process X and lower boundary A. 

Theorem 3.1. Let {Xt} be a Levy process and the loss rate defined in (0). // 
J-^ yv{dy) = 00, then l^ — 00, and otherwise 



l^ = EATi 



1 



K - EAo Jo 



1 



xTTKydx) — 



{K - EAo) 



2{K - EAa) 



(15) 



2{K - EAo) Jo 



K 



(fKix, y, z)v{dy)TTl^{dx)adz), 



where 



-{x — z)^ — 2y{x — z) if y < —X + z, 
'PK{x,y,z)^{ y^ if ~x + z <y < K - X, 

2y{K ~ x) - [K - x)'^ ify>K-x. 
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Proof. The first claim follows immediately if we note that for yv{dy) = oo and 
L > K we have: 

pK poo /"OO 

F> TTKidx) [y - K + x)v{dy) > {y - K)v{dy) = ^. 

Jo JL JL 

The idea of the proof of the second part of the theorem is based on two steps: 1) 
expanding all terms on the right-hand side of and 2) sending L to infinity and 
then a X 0. 

Step 1. In particular, the first term on r.h.s. of ([T4| equals: 



K 



'i°'^'^(\y\>L)v{dy)TiK{dx) 

(pK poo ^ 
9l~ X li^y\>L)v{dy)TTK{dx) 
Jo J ~oo J 

\9l J XTTK{dx) + — + J —y{dy) 

'^{\y\>L)v{dy)iiK{dx) +o(a^). 



K ^2 /-oo 



Similarly, 



/oo 
yl{\y\>L)i^{dy), 
-oo 



(AL^f^^ r TTKidx) f {y-K + xfu{dy) 



0<s<l 



K-x 



and 



a 



2 pa pK f — x-^z 



0<s<l 



2 k 



{x + y — z)'^ ^{dyjn^ {dx)£^{dz) . 



-L 



For X > denote v{x) = v{{x,oo)) and similarly for x < let v_{x) = oo, a;)). 
Then: 



= — c 



K 



yv{dy) ~ a I xTrfc{dx)i^{—L) + i 



K 



{dx)^{dz)u{-L) 



and 



tK 



a / y^idy) + a K yv{dy) 

J L J L 

+{a + a^K) i XTTK{dx)v{L)- KV{L) 
'o 
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The other terms equal: 



2 7^2 rK poc 

e"^E ^ = {I + aK + ^-)V{L) - e"" / e"V(dy)7rK(d2;) + o(a") 

0<s<l "^O "^^ 

and 

E ^ e"^-^(l-e-"^^^) - E ^ (1 - 6""^^^ ) + a^E ^ ^ AL^ + o( 

0<s<l 0<s<l 0<s<l 

with 

^ = k[-L)- e"" / e''yv{dy)nK[dx) 



a') 



0<s<l 



/ e"^ / e''yiy{dy)TTl^{dx)^{dz) 







2 







Step 2. If we now rearrange all terms of ([T4|) using above identities and let i — > cx) 
(note that 9l EXi as L oo) we get: 

pK 2 2 

a(l - e"-^ + aEAo)?^ = -EXia^ xTTK{dx) ~ '^-^ 





2 /-a ^/^ fK—x 

2 



?;'i/(d2/)7r^(da;)C(rf2;) 

J z J —x+z 



2 /-a ^/f ^oo 

2 



((a; - if)^ + 2y{x - K))i^{dy)TT^K{dx)^idz) 

v/^ JK-x 
2 /-a ^if ^— x+z 

((x - z)2 + 2?;(x - z))iy{dy)Trl^{dx)^{dz) 

s J —oo 



a 

,.2ti? 4 h? \^ I „/'„.2 



WEAoEXi +o(a^). 

The proof is completed by dividing both sides of above equation by a(l — e"^ + qEAq) 
and sending a to 0. 

□ 

Assume now that there exists 7 > (7 G 0) such that K{-y) — 0. Define now the 
new probability measure: 

dF-^ 



for which we have E^ATi = k (7) > since on P''' process X is a Levy process with 
Laplace exponent Kj{a) — K(a + 7). We also need two passage times: 

T^ix) = inf{i > : ATt > A'^ + x}, tZ^ ^ inf{t > : ATt < -z}. 

Further, let t^(x) = inf{t > : ATt > If + x}, where A\ = A'^ ^{dy) and 
B"^{x) = Ar^A(j.) — X. The second main result concerns the asymptotics of as 
K ^ 00. 
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Theorem 3.2. Assume that there exists 7 > (7 e 6j such that k;(7) = and 
AC (7) < 00. Then there exists random variable B-^{oo) such that: 

lim ET'e-T^''^^) = ]ETe-T^''(°°\ (16) 

X— >-CX) 

Furthermore, there exists finite constant D such that: 

~ De-^^, as K ^00, (17) 
where we write f{K) ~ g{K) when lim;f^.oo f{K)/g{K) = 1, and 



D = -EXiC^+]ETe-T^''(°°)y eT^PT(Tr^ = oo)y (l - e^^^"^) j zv(d?/) 



da; 



+ / (2/ + 7"'(l-e^")) Kc^y) 

noo naAx r—X-\-Z . . 

+ y y F{Tf{x)<oo)j (1 - e'^'^^'+y-^^j u{dy) ^{dz) dx (18) 
Proof. The proof is based on the following observation: 

' -e7^-c/^ + ^^^^^^ + ^^^^^^ + ^^^^^^"e«^-C,' ^^^^ 
where 

-^1 = / / {y - K + x)v{dy)-KK{dx), 
Jo Jk-x 

h= I I I {x + y - z)u{dy)'K^K{dx)^{dz), 

Jo Jo J -00 

h= / {l-e'^^y-''+^^)v{dy)T,K{dx), 
Jo Jk-x 

h = r r I '^\l-e'^^-+y-'Mdy)TT'j,{dx)adz). 
Jo Jo J -00 

Indeed, note that from equation (7) taking a = 7 we get: 

= 7C^Z^ - 7e^^/f + C^E ^ (1 _ 6"^^^^) + e^^E (l - e^^^^), 



0<s<l 0<s<l 



where we used fact that E eJ'^'dLf = C^/l^. Let e > 0. Wc split ALf into two parts, 
A'^Lf and A^Lf-, corresponding to AX^ G [0, e] and AXg G (e, 00), respectively, and 
ALf into A^Lf and A^Lf, corresponding to AXg G [-e,0] and AXg G (-00, -e), 
respectively. Now we have: 



2 



e-^^E ^ (l-e-A^-^)=e-^ -7(if-i;^)-^E ^ (A^Lf )^ | + o(a, 

0<s<l V 0<s<l 
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0<s<l 0<s<l 

Thus, 

0<s<l 0<s<l 
0<s<l 0<s<l 

Me'). 

Using the fact that — — EX\ we have: 

Q<s<l 

-e^^E Y (l-e^^'^") + ^C,E ^ {A'Lff + ^E ^ (A'Lff 

0<s<l 0<s<l 0<s<l 

If we send e ^> we get: 

1^{C^ - e'^^)7 = 7C^EXi - jC^h - -ie^^h - C^h - e'^^h- 

Now (dH) follows by dividing by (C^ - e''-^)7. 

Note that Ii and are the same like those in [2] Theorem 3.2] and I2 and have 
just additional integral over ^{dz) and we should take x — z instead of x under the 
integrals sign. Using thus the same arguments like in the proof of ^ Theorem 4.1] will 
now completes the proof when we prove weak convergence (jl6p . To this one can use 
classical renewal arguments applied to the process {Xt-a(£j„), ri £ N} (by considering 
ladder height lines being At starting from invariant measure shifted by a from the 
previous position of ladder process). 

□ 

Example 3.1. For a stable M/M/1 queue, that is for Xt — X^iL'i with {(Ti}{i>i} 
being i.i.d. exponentially r.v.'s with intensity and Nt being a Poisson process with 
intensity A < /u, we have vldx) — fiXe~'^^ dx and P''{tZx — 00) = 1 — e"^"^, where 
7 = — A, since considering X on is equivalent of exchanging intensities of arrival 
and service processes. Moreover, choosing saw-like lower boundary given in ([5]) by lack 
of memory of exponential distribution on P''' we have B^{oo) = ei — Y , where Y has 
uniform distribution S,{dx) — ^ {x & [0,a]) and Ci is exponential r.v. with intensity 
A. This gives: 

a \ J fi ^ 
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